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, The electrical properties of a two-dimensional packing of metallic beads are 

studied. Small mechanical perturbations of the packing leads to giant electrical 
^ I fluctuations. Fluctuations are found to be non-gaussian and seem to belong to 

Levy stable distributions. Anticorrelations have been also found for the sign of 
I . these fluctuations. 
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The granular state of matter exhibits fascinating physical properties such as 
memory effects [Q , phase segregation ||^ , heterogeneous mechanical stresses Q , 
etc... Although many experiments have been performed for studying mechanical 
and geometrical aspects of the granular matter Q , only very few reports |Q, ||, ^ 
can be found in the scientific literature for describing electrical properties of such 
systems. In 1890, Branly reported Q that the electrical resistance of packed 
metallic grains is affected by an incident electromagnetic wave. The Branly's 
coherer was the first reported antenna. He discovered also that a small vibration 
of the packing affects the electrical resistance. More recently, Ammi et al. 
measured a non-Hertzian behavior of the electrical conductivity when a packing 
is submitted to high uniaxial pressures. In an other experiment, Vandembrouck 
et al. visualized the electrical paths in a packing of metallic beads for different 
injected current densities. 

In the present letter, we report electrical measurements of a 2d packing of 
metallic beads submitted to small perturbations. Electrical fluctuations are 
statistically analyzed. A physical interpretation is given within the framework 
of Hertzian contacts. Comparison with recent experiments and models for stress 
fluctuations are also emphasized. 

Our experimental set-up is illustrated in Figure 1. The main part of the 
system is a tilted plane of epoxy of 25 x 30 crri^ on which P6O2 beads have 
been placed. A small tilt angle (typically 9 — 10°) can be adjusted by two 
screws. This plane has been made using printed circuit technique allowing to 
design easily the configuration of the electrical contacts. In our experiments, 
three contacts (A, B, C) have been used as described in Figure 1. Contacts 
are placed on specific beads. A plexiglas protection plane has been placed just 
above the epoxy plane to ensure the quasi-2d characteristic of the system. The 
mean diameter of the beads is 2.35 mm and a polydispersity of 2 % has been 
measured. About 3000 beads have been placed on the epoxy plane. Above the 
protection plane, a wooden hammer driven by the parallel port of the computer 
can hit the epoxy plane at the bottom of the packing. 

The experiments have been realized in a constant current regime. A current 
source injects a current i through A and C. The tension U between the point 
B and C is then sent to the computer by a nanovoltmeter. This method is the 
so-called "3-points measurement" . In each series of measurements, the following 
procedure has been repeated a large number n of times: (i) the beads are placed 
in between the epoxy and plexy planes, (ii) the hammer hits the bottom of the 
packing every 10 seconds, (iii) the computer waits 3 seconds and then stores 
the tension between B and C at sampling rate of 0.2 seconds during the next 7 
seconds. 

Figure 2 presents a typical evolution of the voltage [/ as a function of the 
shock number n. The inset illustrates the tension at the scale of a few seconds: 
shocks events are denoted by arrows at that scale. When a shock occurs, U 
is seen to be drastically affected and thereafter sticks at that new value. One 
observes also that there is apparently no systematic behavior at each shock: (i) 
the voltage fluctuation is either positif or negatif, and (ii) the jumps/drops of U 
seem to be characterized by a broad distribution of amplitudes. A compaction 
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of the packing is also observed at the early stages of the experiment, and there- 
fore one observes that the voltage U globally decreases since more and more 
beads are touching. After n = 2000 shocks, we have observed no macroscopic 
drift of the electrical voltage in our measurements! After t = 10000 shocks, 
the packing is still disordered. Indeed, numerous departures (defects) from an 
ordered structures are observed. 

Fluctuations AC/„ are defined as differences of the electrical voltage before 
Un and after Un+i the shocks. Since i is constant during the experiment, fluc- 
tuations AU are mainly due to changes of the conductivity of the packing, i.e. 

Act 

AU = iAR!v-i— (1) 

from the Ohm's law U — Ri. From Eq.(l), one expects that the mean size of 
the fluctuations (AU) would be proportional to the injected current i. In order 
to compare our different data series, we normalized AU by i. Figure 3 presents 
the distribution of fluctuations P{AR) in a semi-log plot for different values 
of the injected current i. So-called "fat tails" for rare events are observed in 
Figure 3, the ocurrence of large fluctuations is several orders of magnitude more 
frequent than those for gaussian tails. The fluctuation histograms arc however 
well described by a Levy stable distribution which has the following behaviors: 

P(|Ai?|) ~exp(-7|Ai?|«) for |Ai?| « I/7 (2) 

and 

Pi\AR\) - |Ai?r("+i) for \AR.\ » I/7 (3) 

where < a < 2 is the so-called Levy index and 7 > is a quantity related to 
the width of the distribution (for a > 1). When a reaches 2, the distribution 
becomes a Gaussian distribution. When a. = 1, the Levy distribution reduces 
to a Cauchy distribution. When a < I, the second moment of the distribution 
diverges. The continuous curves in Figure 3 are fits with Eq.(2) giving a value 
a = 0.9 ±0.1 for both cases. Figure 4 presents in a log-log plot the fat tail of 
the distribution P{\AR\). Power law behaviors (Eq.(3)) are observed. We have 
found a = 0.95 ± 0.10 in good agreement with the a value of the first fits using 
Eq.(2). Both fits in Figures 3 and 4 confirm the Levy scenario for electrical 
fiuctuations. This value is quite low (lower than a Cauchy) and expresses that 
the voltage evolves erratically with frequent large drop and burst excursions. It 
is also worthwile to point out that the Levy index a seem to be independent of 
the injected current i ranging from 3mA up to 100mA in our different series of 
measurements. 

Let us also investigate the temporal correlations between two successive 
fiuctuations. Figure 5 presents the phase space (A[/„, A[/„_|_i) for a current 
i=50mA. The data points are mainly dispersed in an ellipsoid shape which 
is tilted with respect to the horizontal and vertical axis. Moreover, branches 
are observed along the vertical, horizontal and diagonal directions. The latter 
branches indicate the existence of strong correlations for large fluctuations. The 
continuous line is a linear fit through the entire cloud of data points. The slope 
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a = —0.45 ±0.01 of the linear regression Af7„+i = alSUn + b is significantly neg- 
ative, and expresses that two successive fluctuations are highly anticorrelated. 
The slope a is found to be weakly depend on i. 

Finally, it should be noted that Ni beads which are more conducting than 
P6O2 beads lead to similar results: Levy-like distributions and anticorrelations. 
We checked also that the results are not sensitive to the tilt angle 9. 

If the properties of the fluctuations characterizing the phenomenon are found 
to obey non-trivial laws, they should certainly arise from fundamental phenom- 
ena. A physical interpretation of these giant and non-Gaussian fluctuations 
should be found in the framework of the geometry of electrical paths. Indeed, 
the existence of multiple tortuous conducting paths have been reported in 
We conjecture that the fluctuations observed herein are most probably related 
to the fluctuations of stress chains in the packing. Also, the strong anticorre- 
lations that we put into evidence should be related to memory effects in the 
packing. 

As recently underlined by deGennes 0| , both macroscopic electrical conduci- 
tivity a and macroscopic elastical modulus /i of the packing are closely related 
if Hertzian contacts are assumed. This analogy between mechanical strain and 
electrical conductance implies that the distribution of the fluctuations measured 
herein could be related to the stress fluctuations in the packing. Giant stress 
fluctuations have been recently reported in compression Q and photoelastic ex- 
periments 1^ . Cooperative rearrangements such as nucleation of cracks and the 
propagation of these cracks in the packing can lead to stress flcutuation of high 
amplitude. It has been put into evidence that theses cracks have a transient life 
and heal after the displacement of blocks or single beads From the simu- 
lation point of view. Coppersmith et al. ||l^ introduced the so-called g-model 
deducing a probability distribution of vertical forces varying as /exp (—//(/)). 
Howerer, the g-model and more elaborated variants lead to Gaussian distribu- 
tions of stress fluctuations. At our knowledge, the only one model which is able 
to produce rearrangements of force paths at all scales is the so-called Scalar 
Arching Model (SAM) introduced by Claudin et al. 0. The SAM gives rise 
to a broad distribution of the apparent weight variations ISWa measured at 
the bottom of a silo. Claudin has found ||l^ a universal power law behavior 
P(|Aiya|) ~ |AWa| whatever the values of solid friction coefficients. This 
power law behavior is consistent with our results but the relationship between 
our exponent a + 1 and the Claudin's exponent 0.94 is non-obvious. 

In summary, the careful measurements of the electrical properties of a quasi- 
bidimensional packing of P6O2 beads have put into evidence the non-gaussian 
character of conductivity fluctuations. Those non-Gaussian electrical fluctua- 
tions are related to stress fluctuations in the packing. Our measurements confirm 
some predictions of the Scalar Arching Model. 
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Figure Captions 

Figure 1 - A schematic illustration of the experimental setup. Contacts A, B 
and C on beads. Epoxy (E) and Plexiglass (P) Planes which are separated by 
the bead diameter and which are tilted by an angle 9. 

Figure 2 - The electrical fluctuations of U{n) in the quasi-2d packing of P6O2 
beads. The inset presents the evolution of U at the scale of a few seconds. 

Figure 3 - Semi-log plot of fluctuation distributions P{AR) for two different 
current intensities i. The continuous curves are fitted Levy-like distributions 
using Eq.(2). 

Figure 4 - Log-log plot for the tail of P(|Ai?|). The lines are power law fits 
using Eq.(3). 

Figure 5 - Phase space {AUn, AJ7„_|_i). The continuous line is a linear regres- 
sion A{7„_|_i = aAUn + b. 
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